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The  reverse-thrust  performance  of  two-dimensional  and  a\is>mmetric  tuned  deflectors,  ssmmetricallv  situated 
with  respect  lo  the  impinging  jet,  is  analyzed  through  the  use  of  l^si-f  isita's  method,  Belolserkossky's  integral 
method,  and  the  finite  element  method,  lire  predicted  results  are  shown  to  be  in  reasonable  agreement  with  the 
available  experimental  data.  The  real-fluid  effects  such  as  Coanda  effect,  entrainment,  and  the  pressure  ratio 
are  discussed  in  the  light  of  the  present  imiscid  How  analsses. 


Nomenclature 

a,  =«  coefficients 

•I,  » areas  of  suhtriangles  in  a triangle 

.4’  - area  of  a triangle 

J = jet  width  or  diameter 

</{r)  « reverser  geometry 

It  m length  of  the  normal.  Fig.  6 
lilol  = a function 

i - I - II1  1 

/((.'ll  m functional  of  </> 

k = a parameter 

ii  - normal  axis 

p = pressure 

q m magnitude  of  total  velocity 

r * radial  distance 

r,  = nozzle  radius 

R » radius  of  curvature 

R.  = reverser  radius 

s = axis  along  the  free  surface 

,v„  * nozzle  setback  distance 

SA,jm  - element  matrix  for  element  m 

SLA,m  = load  matrix  for  element  w 

i = a parameter.  Fig  I 

ii.  i = velocity  components 

I = free  streamline  velocity 

l„  = normalized  uniform  velocity  in  the  nozzle 

It'  = complex  potential 

v.y.r  = coordinate  axes 

/(  = deflection  angle 

J « slope  of  the  reverser  boundary.  Fig.  6b 

; = a variable 

”,  = area  coordinates 

)/„  = reverse-thrust  ratio 

tl  = direction  of  velocity  vector 

/ = a multiplier 

l>  = density  of  fluid 

it  = a variable.  Fig,  I 

ip  = potential  function 

ip  = stream  function 

i t - a complex  function 


Introduction 

THE  fluid  dynamics  of  a free  air  or  liquid  jet  impinging 
on  a rigid  or  deformable  surface  poses  problems  ol  special 
interest  in  the  analysis  of  impulse  machinery,  thrust  reversers, 
flip  buckets  on  spillways,  rock  and  metal  cutting  with  high- 
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speed  water  jets,  false-tw  isting  of  yarns,  etc.  A complete  discussion 
of  the  mechanics  of  the  flow  situation  would  include  the 
compressibility  and  gravitational  effects,  boundary-layer  effects 
at  the  solid  surface,  and  the  entrainment  effects  at  the  free 
surface,  but  in  most  cases  these  are  relatively  unimportant,  and 
it  is  sufficiently  realistic  to  assume  the  fluid  to  be  inviscid  and 
incompressible.  Even  with  this  simplification,  the  prediction  of 
the  characteristics  of  the  deflected  jet  may  still  pose  exceedingly 
complex  problems  particularly  for  jets  impinging  obliquely 
upon  prescribed  boundaries.  For  example,  the  oblique  impact 
of  a round  jet  on  a plane  surface  has  not  yet  been  analyzed.1 
Thus,  the  analysis  is  often  restricted  to  two-dimensional  or 
axisymmetric  situations  and  recourse  is  usually  made  to  approxi- 
mate methods  of  analysis  where  a free  surface  is  assumed  and 
its  suitability  is  assessed  from  approximate  potential  solutions 
obtained  through  the  use  of  various  numerical  techniques. 

The  two-dimensional  counterpart  of  the  jet-deflection  problem 
has  been  treated  through  the  use  of  the  powerful  analytic- 
function  theory  and  successive  conformal  transformations. 
Sarpkaya2  solved  the  U-shaped,  two-segment,  deflector  problem 
where  the  turning  angle  between  the  segments  is  limited  to  lK)  ‘. 
Tinney  ct  al.3  extended  this  analysis  to  the  rase  where  the 
turning  angle  between  the  symmetrically  situated  segments  is 
greater  than  90  . Later.  Chang  and  Conly4  presented  an  analysis 
for  a bucket  composed  of  a series  of  segments  of  arbitrary 
number,  lengths,  and  angles.  However,  the  basic  as  well  as 
practical  problem  of  the  direct  and  exact  analysis  of  jet  deflection 
by  curved  buckets  remained  unsolved. 

The  three-dimensional  counterpart  of  the  jet-deflection 
problem  has  not  been  solved  in  any  generality.  Attempts  to 
formulate  an  exact  solution  have  been  mostly  unsuccessful  even 
for  axisymmetric  inviscid  flows  with  no  body  forces.  The  case 
of  a circular  jet  striking  a plate  normally  was  analyzed  by 
Schach5  using  approximate  methods  similar  to  those  of  Trefft/1’ 
with  successive  adjustment  of  the  free  streamlines.  Jeppson7 
applied  the  finite  difference  technique  to  the  solution  of  two, 
axisymmetric,  potential-flow  problems,  namely  to  that  of  flow 
from  a nozzle  and  the  cavitating  flow  of  a iet  past  a body 
of  revolution.  Other  noteworthy  contributions  to  the  analysis 
of  the  jet  efflux  from  nozzles  and  orifices  were  made  by 
Southwell  and  Vaiscy.'1  Rouse  and  Abul-Fctouh.1'  Garabcdian.10 
and  Hunt"  through  the  use  of  the  relaxation  and  finite- 
difference  methods.  Schnurr  ct  al.12  used  the  relaxation  method 
to  analyze  the  turning  of  two-dimensional  and  axisymmetric 
jets  from  curved  surfaces  where  the  jet  was  assumed  to  leave 
the  deflector  exactly  parallel  to  the  tangent  at  the  lip  of  the 
deflector  surface.  The  consequence  of  the  difference  between  the 
actual  deflection  angle  and  the  said  tangent  to  the  deflector 
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Fig.  I Transformation  planes;  a)  physical  plane,  b)  If -plane,  r)  /-plane, 
and  d)  C -plane. 

mentally.  Thus,  their  analysis  docs  not  constitute  a solution  to 
the  thrust-reverser  problem 

Ryhming1 3 studied  the  problem  of  a steady  liquid  jet  impacting 
on  a plane  and  on  an  infinite  axisymmetric  cavity  by  employing 
Belotscrkovsky's14  integral  method  As  in  the  case  of  Schnurr  et 
al.,  the  results  of  this  work  do  not  apply  to  the  thrust-reverser 
problem  where  there  are  two  free  surfaces  and  where  the 
asymptotic  slope  of  the  deflected  jet  is  not  parallel  to  the 
tangent  at  the  lip  of  the  deflector  surface. 

In  the  present  work  three  methods  arc  employed  to  investigate 
the  deflection  of  inviscid,  incompressible  curved  surfaces:  Levi- 
Civita's  method.  Belotscrkovsky’s  integral  method,  and  the  finite 
element  method.  The  first  method  is  applied  to  plane  jets 
impinging  symmetric-ally  on  two-dimensional  curved  surfaces 
whose  shapes  arc  specified  in  terms  of  a given  jet-deflection 
angle  and  the  angle  of  departure  at  the  lip  of  the  two- 
dimensional  bucket.  The  other  two  methods  are  applied  to 
straight  circular  jets  impinging  on  axisymmetric  curved  surfaces. 

Plane  Jets  and  Two-Dimensional  Curved  Deflectors 

The  past  as  well  as  the  present  analyses  have  emphatically 
demonstrated  that  almost  all  of  the  difficulties  of  the  analysis  of 
jet  deflection  from  curved  obstacles  are  ascribable  to  the  deter- 
mination of  the  initially  unknown  free  streamlines.  In  view  of 
this  fact  a direct  analysis  of  the  problem,  that  is.  the  determina- 
tion of  the  jet  deflection  angle  and  the  position  of  the  free 
streamlines  via  either  purely  analytical  or  partly  analytical  and 
partly  numerical  techniques  is  at  best  a difficult  problem. 
As  propounded  by  BirkhofT  and  Zarantonello,15  what  one  can 
hope  for  is  an  indirect  solution  of  the  problem,  i.e  obtaining  a 
family  of  obstacle  shapes  which  will  yield  the  prescribed  jet- 
dctlcction  angle  and  provide  a familiarity  between  the  shapes 
so  obtained  and  the  range  and  change  of  direction  of  the  free 


parameters  involved  in  the  analysis.  This  holds  true  for  both 
the  plane  and  axisymmetric  flows.  Should  one  adopt  a purely 
numerical  procedure  (e.g.  finite  difference,  relaxation,  finite 
element,  markcr-and-cell  technique),  it  is  then  and  only  then 
possible  to  begin  with  a prescribed  obstacle  shape  and  to 
approach  a unique  solution  through  successive  iterations. 

It  is  well  known  from  the  free-streamline  theory  that  whenever 
the  solid  boundaries  are  composed  of  straight  segments,  then 
the  ft-plane  also  is  composed  of  straight  lines  (sec,  e g..  Ref  16). 
Since  such  a polygonal  boundary  can  always  be  transformed 
onto  a t-planc  through  the  use  of  Schwa  r/-ChristofTel  transforma- 
tion, all  two-dimensional  jet-deflection  problems  of  this  nature 
may  be,  at  least  theoretically,  solved.  However,  whenever  both 
the  magnitude  and  the  direction  of  the  velocity  vary  along  a 
rigid  boundary,  then  the  fl-plane  is  not  in  general  composed 
of  straight  lines  or  of  lines  intersecting  each  other  at  suitable 
angles  which  would  yield  intcgrable  transformations.  Con- 
sequently, one  will  cither  seek  other  methods  of  handling  the 
curved  boundaries  or  completely  abandon  the  direct  approach 
of  obtaining  a solution  for  a given  geometry.  The  indirect 
approach  seeks  a family  of  solutions  for  a set  of  jet-deflection 
and  departure-angle  conditions  and  then  lets  the  designer 
choose  the  one  best  suited  to  his  needs. 

A relatively  direct  approach  consists  of  rounding  ofT  of  the 
corners  of  buckets  otherwise  comjsosed  of  straight  segments. 
Generally,  this  procedure  results  in  a considerable  complication 
of  the  mathematical  problem  and  requires  suitable  analytical 
and  physical  simplifications.  The  existing  methods  may  be 
classified  as  follows.  Rounding  olT of  the  corners  a)  by  assuming 
a constant  pressure  transition,17  b)  by  assigning  a special 
pressure  distribution."'  c)  by  modifying  the  /-plane  and  the 
corresponding  terms  in  the  Schwarz-Christoffel  transformation ,g 
(a  brief  discussion  of  this  method  and  the  difficulties  associated 
with  it  are  also  described  by  Carrier  ct  al.).30  and  finally, 
d)  through  the  use  of  special  transformations  such  as  Riemann- 
Hilbcrt  transformation.31-23  hodograph-plane  transforma- 
tions.33 etc.  As  noted  previously,  some  of  these  methods  require 
the  solution  of  complicated  differential  or  integral-equation 
problems  and  some  consist  in  first  solving  the  problem  with 
sharp  corners  and  subsequently  replacing  the  original  boundary 
with  a modified  one  having  rounded  corners.  Often  the  time 
and  effort  required  to  round  off  the  sharp  corners  may  not  be 
commensurate  with  the  need  for  a solution  and  it  may  be 
desirable  to  adopt  an  indirect  approach  where  the  shape  of  the 
curved  boundary  is  not  initially  known  but  obtained  as  part  of 
the  solution. 

l>evi-Civit»'s  Method 

Lcvi-Civita1  5-24  may  be  said  to  have  solved  the  inverse 
problem  of  describing  the  class  of  all  two-dimensional  jets 
divided  by  curved  barriers.  Although  no  specific  solutions 
have  been  presented,  the  problems  of  determining  such  flows 
were  reduced  in  special  cases35-36  to  the  solution  of  nonlinear 
integral  equations  with  appropriate  boundary  conditions.  The 
method  has  not  previously  been  applied  to  the  solution  of  the 
deflection  of  a two-dimensional  jet  from  a finite  two-dimensional 
curved  boundary.  It  consists  of  the  definition  of  a complex 
function  id  as 

id  = O + iLniq/l))  (I) 

and  the  definition  of  the  shape  or  the  curvature  of  the  obstacle 
in  terms  of  the  coefficients  of  a polynomial  representing  id. 

For  the  problem  under  consideration,  the  physical  c-plane 
(see  Fig.  la)  consists  of  a jet  of  width  2 <1  and  velocity  \) 
impinging  upon  a curved,  symmetric,  two-dimensional  bucket 
BAB'  The  complex-potential  function  IV'  is  given,  as  usual,  by 
U'  = ij>  + iip.  A straightforward  application  of  the  Schwar/- 
Christoffcl  transformation  establishes  a relation  between  the  r- 
and  H'-plancs.  Thus,  one  has 

U'  = — (lj//7!)[/ai(l— r)+/./i(l+/)J  (2) 

The  parameter  k in  the  /-plane  is  uniquely  related  to  the 
deflection  angle  It. 
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A typical  representation  of  0(a). 


l:>(  -t  I ) = — 71.  (ltd)  ^ ±71 '2  ■ 100,  I «(  t)  ss-  71,  ll>(0)  — () 

The  relationship  between  0{n)  and  a may  be  chosen  at  will. 
Let  this  relationship  be  repre-ented  by  a simple,  symmetric 
variation  as  shown  in  Fig.  2.  Equation  (10)  may  now  be 


written  as. 

COSIT 

itj(rj)  = Hji  + ih(a)  + 0(n) 

(12) 

1 -t  sin  it 

where 

(i(it)  = £u„sinnn,  n=l.3, 5. ... 

(13) 

and 

0(a)  = -(ir/2)+  Y. a* cosn  it 

(14) 

The  coefficients  a,  are  giun  by 


The  upper  half  of  the  f-plane  is  now  transformed  into  the 
inner  region  of  a semicircle  of  unit  radius  (see  Fig  Id)  through 
the  transformation  given  by  Lcvi-Civita24  as 

r '12  <3> 

The  circular  arc  in  Fig.  Id  represents  the  solid  boundary 
over  which  the  direction  of  the  velocity  vector  is  assumed  to 
be  prescribed  and  the  line  BCDCB’  represents  the  free  steam- 
lines where  the  magnitude  of  the  velocity  is  known. 

The  functions  IF  and  «>  may  be  written  as 

-(|/fj)|t/H'/i/--)  = (q  l')e  * (4) 

and 


/ 1 </IF\ 

.v  i r 

(5) 

d n.L 

which,  in  terms  of  reduces  to 


a„  = [2/(it,  - ir,|n2](cosn  it,  - cosn  ir2)  (15) 

in  which  it,  and  n2  represent  the  two  parameters  shown  in 
Fig.  2.  l or  m = 0,  a0  = - n/2  as  expected  from  the  average  value 
of  the  function  0(a). 

To  complete  the  analysis  one  needs  to  evaluate  the  term 
exp  [iVu(rr)]  which  appears  in  Eqs.  (7)  and  (8).  Writing, 

t, = [,  | + sin  cos  ,,  ( 1 6) 

and  separating  the  real  and  imaginary  parts,  one  has 

S'*"'  = (I  •+■  sin  irFcosire  '"'"[cos  (((it)  + i sin  0(a)]  (17) 

Combining  Eqs.  (7).  (13),  (15k  and  (17)  one  has 

.v  1 f"  ,2k2(sin  it)(I  +sin  a)  . 

= e , sin  (Ho)do  (18) 

d n\,  2 I -k  cos2  a 

and 


(6) 


Defining  ; = re\p(i>r)  and  noting  that  over  the  bucket  only 
the  argument  and  along  the  free  streamlines  only  the  modulus 
varies,  the  combination  of  Eqs.  (2)  and  (6)  yields. 


and 


d 


d 


„ 1 |>.,2*,‘^*»*inrT 

7i  I 1 - A;2  cos2  a 

I f , , , k2(r-r  2)/2r 
rcj  k2(r  + r V/4- 


da 


. dr 


(7) 

(8) 


where  Eq.  (7)  gives  the  coordinates  of  the  bucket  and  Fq.  (8) 
the  coordinates  of  the  free  streamlines 
The  function  w is  assumed,  according  to  Levi-Civita,  to  be  of 
the  form 


oil-,  I — //-«[(/ — s)/(i  + + </,,+  (/(s  4-  it2.,2  + • (9) 

where  the  logarithmic  term  can  be  identified  as  that  correspond- 
ing to  a free  jet  impinging  normally  upon  a flat  plate.25  The 
additional  terms  i.i  the  polynomial  whose  coefficients  are  to  lie 
prescribed  later  modify  the  expression  for  v>  and  give  a suitable 
camber  to  the  rigid  boundary. 

Noting  that  u2  <i4  <;h must  be  taken  equal  to  zero  because 

of  symmetry  of  the  (low  and  that  r = I over  the  solid  boundary 
and  it  = 0 or  a = n over  the  free  streamlines.  Eq.  (9)  may  be 
reduced  over  the  bucket  to 

Icosit'  n 

w(ff)  = i/jt  ~ 7 - ± „ +<(,  vl"  + a}ei,"  + iin->‘"+  ...  (10) 

I + sin  it  2 


where  plus  sign  is  to  be  used  for  n > a > it/ 2,  and  minus  sign 
for  it  2 > it  > 0.  Over  the  free  streamlines,  Eq.  (9)  reduces  to 

2 r 

f.o(r)  = —tan  , +</i r + tiKr3  + ii;rf+...  (II) 

I - r4 


The  coefficients  nlf  a}. . may  now  be  chosen  to  obtain  different 
bucket  shapes.  In  this  study  a Fourier-series  representation  was 
used  primarily  because  of  the  versatility  afforded  by  it  in 
representing  different  bucket  shapes.  The  evaluation  of  the 
coefficients  requires  a statement  of  the  conditions  at  the  tip  and 
the  axis  of  the  bucket  Let  the  jet  departure  angle  at  the  tip  of 
the  bucket  be  equal  to  ±ir,  and  the  curvature  at  the  stagnation 
point  be  continuous.  Then  one  has 


j-'l\ 

d n\,2 


,2k2(sini7)(l  +sino) 

— r-  cos  0(a)  da  (19) 

I -k  cos  IT 


Equations  (18)  and  (19)  yield  the  coordinates  of  the  curved 
boundary  through  the  use  of  the  appropriate  values  of  it, 
and  a2.  This  in  essence  completes  the  application  of  the  Lcvi- 
Civita’s  method  to  the  analysis  of  the  jet  deflection  from  a 
two-dimensional  curved  bucket  similar  to  that  shown  in  Pig.  1. 

The  deflection  angle  is  evaluated  by  replacing  in  Eq.  (II) 
the  corresponding  value  of  r at  the  point  C i.c., 

r = [I  -(I  -A2)1  2]/A  (20) 

The  foregoing  equations  have  been  numerically  integrated  for 
a series  of  representative  values  of  it,  and  n2.  The  resulting 
family  of  bucket  shapes  are  shown  in  Fig.  3.  A special  case  for 
it,  = 0.31416  and  a2  = 0.47124.  and  k = 0.946.  (/(  = 68  k is 
show  n in  Fig.  4. 

Evidently.  Lcvi-Civita's  method,  with  the  novelties  introduced 
into  it  in  this  study,  is  sufficient  not  only  to  generate  a family 
of  suitable  two-dimensional  curved  reversers  but  also  to  round 
off  the  corners  of  buckets  composed  of  straight  segments.  For 
example,  by  letting  (it,  ~a2)  -*0,  one  finds  that  the  bucket 
shape  generated  by  this  method  exactly  approaches  that 
previously  solver1  by  Sarpkaya.2  For  slightly  larger  values  of 
(it i a2).  the  corner  PQ  (see  Fig.  4)  is  rounded  off.  This  enables 
one.  for  example,  !o  evaluate  the  sensitivity  of  the  deflection 
angle  to  the  rounding  off  of  the  sharp  corner  on  a U-shaped 
bucket. 


Fig.  3 Family  of  two-dimensional  curved  buckets  for  it,  = 0JI4I6. 
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Fie.  4 A ttpical  curu-d  bucket  for  a,  0.31416.  n2  -0.47124,  k 
0.946.  n d - 3.05.  h d » 1.07.  iind  /(  = 6R  . 

The  present  work  may  be  easily  extended  to  the  analysis  of 
the  eases  where  the  lip  angle  of  the  bucket  is  other  than  n 
or  the  bucket  has  a cusp  on  its  axis  of  symmetry  or  to  the 
cases  where  the  impinging  jet  is  not  entirely  free  and  emerges 
from  a prescribed  nozzle. 

Analysis  of  Axisymmctric  Curved  DeHcctors 

As  noted  earlier,  the  three-dimensional  counterpart  of  the 
jet-delicction  problem  has  not  yet  been  solved  in  any  generality. 
The  nnite-dilTerence  and  relaxation  techniques  suffer  in  general 
from  convergence  and  accuracy  problems  and  nearly  all  resort 
to  simple  trial-and-error  procedures  to  locate  the  free  surfaces 
and  to  satisfy  the  boundary  condition  that  the  free  surfaces 
be  stream  surfaces  of  constant  velocity  In  the  following  two 
fairly  new  methods,  namely  Bclotserkovsky’s14  integral  method 
and  the  finite  clement  method  are  used  to  analyze  representative 
axisymmctric  cases. 

Application  of  the  Bclotserko*skj's  Method 

Firstly,  the  coordinates  ,s  and  u along  and  normal  to  the 
unknown  axisymmctric  free  surface  arc  chosen  and  the  flow- 
region  is  divided  into  three  zones  as  shown  in  Fig.  5.  Zone  I 
represents  the  nozzle  and  the  oncoming  jet.  Zone  II  represents 
the  iet  sheet  formed  over  most  of  the  bucket,  and  finally. 
Zone  III  represents  the  deflected  jet.  Zones  I and  II,  and  II  and 
III  are  separated  by  the  stagnation  normal  and  the  tip  normal 
respectively. 

Secondly,  the  functional  forms  of  all  the  integrands  appearing 
in  the  equations  of  motion  are  assumed  and  the  equations  are 
integrated  to  yield,  through  successive  iterations,  the  relation- 


Kig.  5 a is  pica!  axissinmetric  bucket  and  the  definition  of  three 
/ones  of  How  (S  V ; stagnation  normal,  MF  lip  normal). 


Fig.  6 Coordinate  axes  and  geometrical  relations  a)  in  /.one  I, 
b)  in  Zones  II  and  III. 

ships  between  the  geometrical  quantities  describing  the  jet 
boundary  and  the  values  of  the  flow  variables  on  the  boundaries. 

Dividing  all  velocities  by  the  frecstrcam  velocity  Vj.  all 
pressures  by  0.5/iU/,  and  all  lengths  by  the  nozzle  radius  r„.  one 
writes  in  s-n  coordinates  (see  Fig.  6a)  with  corresponding 
velocity  components  ii  and  r the  continuity  and  momentum 
equations  for  a steady,  incompressible,  axisymmctric.  inviscid 
(low  as. 1 •’ 

it  f*  , dll  ( li\ 

^ (rti)</ii  — (ml*  ^ + ( I + - j (rr)*  = 0 (21) 

and 

|-l*-a~ru"  + [(l  + J)rt2r1  + r)Ji 

I (2nJ + p)r//i  + cos(.>  I ( I + 7. ) ptln  = 0 (22) 

J » K Jo'  n / 

where  index  h means  evaluation  of  the  quantity  within  the 
parentheses  or  brackets  at  n = h and  R represents  the  radius 
of  curvature  in  the  meridional  plane  of  the  free  jet  surface. 
In  the  derivation  of  the  above  equations  implicit  is  the  condition 
that  p = r = 0 at  n = 0. 

A first  approximation  to  the  problem  is  usually  obtained 
by  assuming  the  velocities  vary  linearly  with  u,  i.c.. 

n = I -(ii  /i|(l  ii*),  i =(u//i)r*  (23) 

where  ii*  and  r*  are  the  unknown  values  of  u and  r along 
the  axis  of  symmetry  in  Zone  I and  along  the  rigid  boundary 
in  Zone  II  (see  Fig.  6b).  Then  Fq.  ( 2 1 ) at  the  nozzle  exit  reduces  to 

»h  = 3U„-2  (24) 

For  nozzles  for  which  there  is  no  back-pressuring  effect  due 
to  the  bucket  proximity;  i.c  for  large  values  of  s„  (see  Fig.  5). 
one  has  \'„  = u*  = I. 

In  Zone  I.  Eq.  (21 ) reduces  to 

ra/i(2  + ii*)  = 3 (25) 

and  Eq.  (22)  takes  the  form 

(rur)  dn-  f r(2iiJ  + p)i/ii  + cosff  I I I + Ipi/n  = 0 

K (»  J.l  ' Tk 

(26) 
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wind)  reduces  through  the  use  of  ihc  gcomelrieal  relationships 
and  Eqs  (23)  and  (25)  to 

dO 

A , + H = 0 (27) 

(Is 

in  which  -I  and  II  are  functions  of  h,  it,,  and  0 and  were 
previously  given  bv  Kshming  " 

In  Zone  II  I qs.  (21 ) and  (22i  reduce  to 

k^r.+  rj  > M»/itr11+  2z*l  = 3 (2S) 

and 

JD 

r . + q = o (2*>) 

(Is 

in  which  /’  and  Q are  functions  of  h.  u*.  rh,  <i.  and  (I." 

In  Zone  III.  the  equations  derived  for  Zone  II  hold  true 
except  that  p = 0 everywhere,  since  u,  = I and 

II,  = sin  (D-d)  (30) 

r*  = cos  (d-d)  (31) 

Zones  I II  and  III  are  related  by  the  conditions  pertaining 
to  the  stagnation  normal  .S'Y  and  the  tip  normal  FM  (sec 
Fig  5).  Noting  that  at  the  stagnation  point  S.  \ = 0.  rh  = 0. 
and  i//,  = 0.  one  has  from  Fq.  (2X),  for  the  stagnation  normal, 
the  condition 

h2  cos  0 — 1.5  (32) 

Similarly,  noting  that  for  the  case  under  consideration  <5  = n/2, 
(the  tangent  to  the  lip  of  the  bucket  is  taken  parallel  to  the 
v-axis),  one  has,  for  the  tip  normal,  the  condition 

h2(2  cos  0-  cos2  It)  - Virjcos  0-  I ) = 3 (33) 

liquation  (32)  defines  a curve  / 1 = /i(d)  on  w hich  the  end  points 
ofall  possible  stagnation  normals  have  to  lie.  which  passthrough 
the  stagnation  point  S.  Equation  (33)  defines  a similar  curve 
on  which  all  possible  tip  normals  lie. 

The  foregoing  equations  were  used  to  analv/e  the  deflection 
of  axisymmctric  jets  from  buckets  composed  of  a flat  plate 
and  a quarter  circle  as  shown  in  Fig.  5.  The  calculations  were 
performed  as  follows.  First,  a deflection  angle  /( and  the  nozzle 
velocity  l„  were  assumed.  Then  the  thickness  of  the  jet  DE  = </. 
and  the  coordinates  r,  and  r„  were  calculated  from  the  equation 
of  continuity.  Obviously,  d,  r(,  and  r„  arc  subsequently 
recalculated  on  the  basis  of  the  iterated  value  of  1^  The  point 
F.  is  connected  to  the  point  F (the  lip  of  the  bucket)  by  a 
straightlinc  as  a first  approximation  to  the  upper  free  surface 
and  the  line  F.F  is  temporarily  regarded  as  a rigid  boundary. 

The  governing  equations  were  programmed  according  to  a 
standard  Rungc-Kutta  sub-routine  and  then  iterated  as  follows. 
The  calculations  started  at  the  tip  of  the  nozzle  where  0 = 0 
and  it,  = 3F„-  2 and  proceeded  along  the  free  surface.  The 
coordinates  of  the  free  streamline  h and  0 were  calculated  at 
suitable  intervals.  Finally,  a point  was  reached  at  which  Eq. 
(32).  (the  stagnation  normal  condition)  was  satisfied.  Clearly,  the 
stagnation  normal  can  intersect  the  .v-axis  at  a point  other  than 
v=0  and  still  satisfy  Eq.  (32).  Thus,  it  was  necessary  that 
\ = 0 at  rk  = 0.  i.e.  the  stagnation  point  must  be  located  on  the 
bucket.  To  this  end,  the  velocity  V0  was  iterated  upon  until 
both  conditions  [Eq.  (32)  and  \ =0  at  rh  =0]  were  satisfied. 
Physically,  this  iteration  corresponds  to  the  first  order  calculation 
of  the  back-pressuring  effect  of  the  bucket  proximity  on  the 
nozzle  flow  . 

Integration  of  the  Eq.  (29)  in  Zone  II  begins  at  first  with 
the  y„.  ru,  h,  0,  and  R values  calculated  at  the  end  of  the 
Zone  I.  The  free  streamline  calculated  with  these  initial  values 
define  a spline  shift  (see  Fig.  7)  due  to  the  saddle-point  character 
of  the  equations  at  the  singular  point.  Consequently,  a sequence 
of  points  along  the  initial  curve  defined  by  Eq.  (32)  arc  chosen, 
and  the  corresponding  integral  curves  are  calculated.  These, 
in  turn,  determine  an  upper  and  a lower  bound  for  the  proper 
location  of  the  desired  solution.  The  process  is  then  repeated 
until  a sufficiently  accurate  bound  is  determined.  The  above 
procedure  requires  about  ten  iterations.  The  consequence  of  the 
saddle-point  character  of  the  equations  in  the  vicinity  of  the 


stagnation  normal  is  that  there  is  a sector  of  overlap  in  the 
jet,  going  upstream  and  downstream  from  the  stagnation  normal. 
A similar  behavior  has  been  found  previously  with  the 
Helotserkovsky  method  in  similar  situations  using  a first-order 
scheme."  2 Integration  in  Zone  II  is  terminated  at  the  tip 
normal  at  which  Eq.  (33)  is  satisfied. 

The  determination  of  the  shape  and  position  of  the  deflected 
jet  in  Zone  III  requires  three  kinds  of  successive  iterations. 
The  first  is  the  determination  of  the  appropriate  values  of  h and  II 
at  the  tip  normal  which  satisfy  Eq.  (33)  and  yield  a sufficiently 
long  portion  of  the  integral  curve  between  the  upper  and  lower 
bounds.  This  iteration  yields  a lower  free  surface  with  zero 
pressure  and  unit  freestream  velocity  but  the  freestream  condi- 
tions (a,  - l.p  - 0)  at  the  upper  surface  are  not  yet  satisfied 
since  the  upper  surface  was  initially  assumed  to  be  a solid 
boundary  Eurthcrmorc.  the  calculated  stream  surface  does  not 
necessarily  coincide  w ith  the  assumed  lower  free  surface  because 
the  assumed  jet-deflection  angle  is  either  larger  or  smaller  than 
the  actual  deflection  angle.  This,  in  turn,  requires  an  iteration 
on  II.  For  this  purpose,  the  point  I)  is  located  on  the  calculated 
lower  free  surface  and  the  new  values  of  r,  ru.  and  d are  calculated 
from  the  equation  of  continuity.  The  rotation  of  the  upper  free 
surface  constitutes  a first-order  correction  to  the  deflection 
angle  /(.  Then  the  calculations  noted  above  were  repeated  until 
the  deflection  angle  did  not  change  more  than  2 . This  required 
from  five  to  ten  iterations  depending  on  the  initially  assumed 
value  of  the  jet-deflection  angle. 

The  third  kind  of  iteration  was  made  to  determine  the 
position  and  curvature  of  the  upper  free  surface.  As  noted  earlier 
FF.  was  initially  assumed  to  be  a straight  rigid  boundary  and 
thus  the  calculated  tangential  velocities  along  it  were  not  every- 
where equal  to  unity,  particularly  in  the  vicinity  of  the  lip  of 
the  bucket.  Then  the  points  on  the  upper  boundary  were  moved 
inward,  along  the  ii-axis  if  u,  > I and  vice  verse  by  an  amount 

An  = Hu,2-  I)  (34) 

in  which  /.  is  an  assigned  multiplier  and  it,  is  the  previously 
calculated  tangential  velocity.  Experience  has  shown  that  the 
iteration  will  converge  even  for  a very  crudely  assumed  boundary 
for  a value  of  z.  = 0.10.  As  to  the  reason  for  using  In,2  - I),  it 
was  simply  for  the  purpose  of  accentuating  the  error  and 
accelerating  the  correction  of  the  boundary  points  at  which  the 
velocities  most  differ  from  unity.  At  the  end  of  each  iteration, 
the  corrected  boundary  was  smoothened  through  the  use  of  a 
generally  available  CUR  VEFIT  subroutine  by  passing  a smooth 
curve  through  the  calculated  points  and  recalculating  the  inter- 
mediate points.  Then  all  three  types  of  iterations  on  the  lower 
and  upper  boundary  were  repeated  until  the  velocities  everyw  here 
were  within  I +0.05.  The  results  of  a typical  case  calculated  in 
the  manner  described  above  is  shown  in  Fig.  5 in  which  / 1 = 
44  .5  . 


Fig.  7 Spline  shift  in  Zones  II  and  III. 
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Fig.  8 I-  ini  I r element  representation  of  jet  deflection  from  a hemi- 
spherical revrrser. 


It  is  evident  from  the  foregoing  that  the  Belotserkovsky's 
integral  method  has  certain  limitations.  Firstly,  the  assumption 
of  a linear  velocity  along  s and  n coordinates  cannot  be  valid 
for  all  bucket  shapes  and  relative  jet  radii.  The  variation  of  u 
near  the  free  surface  is  rather  large  particularly  for  very  deep 
buckets.  In  such  cases  Belotserkovsky's  method  may  produce 
rather  complex  and  physically  unacceptable  flow  situations  in 
the  vicinity  of  the  stagnation  point.  Although  no  specific  criteria 
can  be  developed  for  the  bucket  shapes  to  which  the 
Belotserkovsky’s  method  can  be  successfully  applied,  it  may  be 
recommended  on  the  basis  of  the  present  study  that  the  bucket 
radius  he  larger  than  1 5 r„  and  the  bucket  depth  be  smaller 
than  2r„.  Clearly,  the  specification  of  the  radius  of  curvature 
alone  at  the  stagnation  point  in  terms  of  the  jet  radius  is 
not  sufficient.  As  to  the  merits  of  the  method,  it  is  rather 
straightforward  and  yields,  within  the  ranges  of  the  parameters 
recommended,  results  which  are  sufficiently  accurate  for  the 
calculation  of  the  pressure  distribution  on  the  bucket  and  the 
reverscr  thrust. 

Application  of  the  Finite  Flement  Method 

Zicnkiewicz  and  Cheung28  proposed  in  1965  the  application 
of  the  finite  clement  method  to  the  solution  of  field  problems 
involving  the  equations  of  Laplace  and  Poisson.  Since  then  a 
significant  number  of  applications  of  the  finite  element  method 
to  fluid  dynamics  has  appeared  in  the  literature. 28  The  method 
has  recently  been  applied  to  several  jet-efflux  problems  involving 
only  one  freestream  surface  and  relatively  small  jet  contraction 
by  Chan30  and  Chan  and  Larock.31 

The  present  analysis  is  devoted  to  a determination  of  the 
angle  of  deflection,  the  location  of  the  freestream  surfaces,  and 
the  velocity  and  pressure  distributions  caused  by  a finite  curved 
bucket  placed  symmetrically  with  respect  to  the  axis  of  an 
axisymmetric,  inviscid  jet  issjing  from  a uniform  nozzle. 

The  governing  equation  is  then  given  by 

<t>.xx  + 4>.Jr  + <t>.r,  = 0 (35) 

For  an  axisymmetric  flow,  the  solution  to  the  Laplace-field 
equation  satisfying  the  specified  normal-velocity  conditions  M.„)* 
is  given  bv  that  admissible  function  <j>  which  minimizes  the 
functional32 

/(4>)  = Jip|”|  [(</>.i)2  + (<A.,)2]rt/r</v  — 2?rp  j"  4>(4>.n)“rils 

(36) 

in  which  A is  the  half  of  a meridional  section  of  the  flow 
and  C is  a portion  of  the  curve  bounding  this  area  where  the 
normal  derivative  is  prescribed.  The  first  integral  in  Eq.  (36) 
represents  the  kinetic  energy  of  the  fluid  within  the  entire 
control  volume  and  the  second  integral  represents  twice  the 
work  done  by  the  impulsive  pressure  ptp  on  the  boundaries  in 
starting  the  motion  from  rest 

The  approximate  minimization  of  the  functional  H4>)  is 
accomplished  by  dividing  the  field  of  interest  into  N triangular 
elements,  with  corner  and  midpoint  nodes,  and  writing 

#(*)  = !/'<*)  (37) 

i 


The  variation  of  </>  within  each  element  is  represented  by  a 
second-order  polynomial  The  use  of  higher  order  polynomials 
leads  to  extremely  tedious  arithmetical  manipulations.  In  terms 
of  the  area  coordinates  = Au'Am  (see  Fig.  8);  this  leads  to 

<t>m  ” C i (2C  i — I +Cj(2^j—  1Mj  + sj(2,j  — IMj  + 

"KiSatta  +4>iv.Ms  + ‘1»iCM(>  (38) 

The  velocity  components  are  given  by 

t\m  = <t\xm  and  i’,"  = <jt,m  (39) 

The  substitution  of  Lqs.  (38)  and  (39)  in  Eq.  (36)  and  the  use 
of  the  Ritz  technique  yield, 

= S V<V  - SLA,m  «*  0 (40) 

in  which  SA,"  and  SLA,m  represent  respectively  the  elemcnt- 
stifTness  matrix  and  the  corresponding  load  matrix  for  a tri- 
angular element  Chan  and  Larock31  who  tabulated  SA^f  and 
SI  A",  assumed,  in  the  derivation  of  the  matrix  SI  A m a 
constant  normal  velocity  on  a boundary  side  of  an  element 
This  condition  which  conflicts  with  the  linear  interpolation  of 
the  velocity  field  has  been  relaxed  along  the  lines  recommended 
by  Kotchergenko  and  Amorim33  and  the  coefficients  SLA,” 
have  been  corrected  accordingly. 

Iteration  Scheme 

As  cited  earlier,  the  deflected  jet  is  characterized  by  two 
initially  unknosv  n axisymmetric  freestream  surfaces.  Thus  special 
attention  must  be  focused  on  finding  a suitable  iterative 
procedure  for  systematically  approaching  the  final  positions 
of  the  boundaries.  The  iteration  is  terminated  when  the  boundary 
conditions  are  satisfied  within  a prescribed  absolute  maximum 
error.  The  boundary  condition  to  be  satisfied  is  that  the  free 
surfaces  be  streamlines  of  constant  velocity. 

The  grid  is  divided  into  two  major  regions.  The  first  region 
consists  of  elements  whose  coordinates  are  fixed  once  for  all 
(see  Figs.  8 and  9).  This  region  is  well  within  the  interior  of  the 
flow  and  is  not  expected  to  be  intersected  by  the  freestream 
surfaces.  The  second  region,  i.c.  the  moving  grid,  consists  of 
elements  whose  coordinates  are  recalculated  each  time  the  free 
surface  is  moved.  The  movement  of  the  grid  is  made  in  such  a 
manner  that  each  element  maintains  a shape  more  or  less 
compatible  to  its  original  shape.  The  grid  inside  the  exiting 
jet  (along  FE  in  Fig.  9)  is  relocated  by  moving  each  middle 
node  to  its  correct  position  between  two  opposite  nodes  on 
each  side  of  the  jet,  (see  Fig.  8).  In  regions  of  high  velocity, 
relatively  small  triangles  and  in  regions  of  rather  low-  velocity, 
larger  triangles  are  used.  The  examp’c  shown  in  Fig.  8 contains 
288  elements  and  683  nodal  points.  The  bandwidth  was  60. 

Iteration  of  the  free  surfaces  began  with  the  assumption  of  a 
deflection  angle  /<  and  the  sketching  of  the  two  free  surfaces 
as  carefully  as  possible  on  the  basis  of  past  experience.  Then 
the  assumed  boundaries  were  regarded  as  rigid  boundaries.  The 
nozzle  length  was  varied  from  2 r„  to  6r„  and  the  jet  length  FE 
from  0.6 r0  to  r„  in  various  programs  with  no  noticeable  difference 
in  the  results.  The  thickness  of  the  jet  DE  = d.  and  the  radial 
coordinates  rt  and  ru  (see  Fig.  9)  were  calculated  from  the 
equation  of  continuity  by  assuming  K = I Obviously,  d.  rt, 
and  rD  are  subsequently  recalculated  on  the  basis  of  the  iterated 
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values  of  The  point  E is  connected  to  the  point  F by  a 
straight  line  as  a first  approximation  to  the  upper  free  surface. 

The  first  run  through  the  computer  calculates  the  velocity  at  C 
(see  Fig.  1).  i.e,  at  the  lip  of  the  nozzle.  It  turns  out  to  he 
larger  than  unity  (the  correct  value  is  \'c  = 1 1,  since  I.  is  assumed 
to  he  equal  to  unity.  K is  immediately  corrected,  to  a first-order 
of  approximation,  by  writing 

F‘=i;‘  »_(f;“  "-1)1'.“  " (41) 

which  is  used  in  all  subsequent  iterations.  Then  r,.  rt>,  and  d 
are  recalculated  as  in  the  case  of  the  application  of  the 
Bclotscrkovsky’s  method. 

The  procedure  which  has  enabled  the  assumed  boundaries  to 
converge  in  a systematic  manner  to  their  final  positions  may 
be  described  as  follows.  Let  the  velocity  at  an  arbitrary  point 
along  the  assumed  boundary  be  I,.  Then  the  boundary  is 
moved  inward,  along  a line  normal  or  nearly  normal  to  the 
boundary,  if  T,  > I and  vice  versa  by  an  amount 

a,i  = ;.(l;j-i)  (4?) 

In  these  calculations,  the  multiplier  z was  taken  equal  to  0.015. 
At  the  end  of  a given  number  of  iterations,  the  velocities 
along  one  or  both  boundaries  may  approach  values  other  than 
unity  because  the  assumed  deflection  angle  is  not  necessarily  the 
correct  one.  Clearly,  the  deflection  angle  adjusts  itself  at  the 
end  of  each  iteration  but  this  adjustment  is  rather  small  (about 
0.2  ).  Thus,  for  a large  correction  in  /(.  say  20',  about  100 
iterations  are  needed.  Instead,  live  iterations  were  carried  out 
for  the  assumed  deflection  angle  and  then  /(  was  incremented 
by  2°.  clockwise  or  counter-clockwise  depending  on  whether 
the  velocities  in  the  upper  free  surface  was  smaller  or  larger 
than  unity  When  the  velocities  everywhere  were  within  I ±0.05, 
/(  was  no  longer  incremented  by  2;  and  the  number  of  iterations 
was  increased  to  25.  The  calculations  were  terminated  when  the 
velocities  everywhere  were  within  1 ±0.015  The  entire  iteration 
for  a given  nozzlc-reverser  geometry  required  approximately 
10  mins  on  an  IBM-360, 67  computer.  The  program  was  written 
in  FORTRAN  and  double-precision  arithmetic  was  used. 

Examples  and  Results 

The  method  and  the  procedures  described  above  have  been 
used  to  analyze  the  characteristics  of  axisymmetric-jet  deflection 
from  hemispherical  (sec  Fig.  8|  as  well  as  fairly  shallow'  (see 
Fig.  5)  target-type  thrust  reversers.  This  type  of  reversers  give 
the  desired  amount  of  reverse  thru  * without  affecting  the  engine 
operation  and  lend  themselves  t stowage  with  a minimum 
amount  of  boattail  or  base  drag. 

The  geometrv  of  the  nozzlc-reverser  combination  is  uniquely 
defined  by  the  ratios  RJr„  and  sjr„  (sec  Figs.  5 and  9)  once 
the  bucket  and  nozzle  shape  are  decided  upon.  It  is  clear,  at 
least  from  the  experiments3'*  that  there  is  a unique  combination 
of  R,ur„  and  sjr„  for  which  the  reverse  thrust  is  maximum. 
The  determination  of  this  combination  and  the  calculation  of  the 
resulting  thrust  constitute  the  essence  of  the  practical  problem. 
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In  the  aircraft  industry,  the  efficiency  of  a thrust  reverser  is 
expressed  in  terms  of  a “reverse-thrust  ratio"  r/„  defined  by  the 
ratio  of  the  actual  reversed -jet  thrust  to  the  forward  jet  thrust 
of  the  nozzle  alone.  In  other  words,  >;*  is  given  by 

>Ik  = (ttpr.2  F„  V)  cos  //)/( nprm*  Vf)  = ( VJVt)  cos  /(  (43) 

Evidently.  * >hi(RJr0.sjrJ  since  both  V,  and  /(  depend  on 
RJr.  and  sjr.  for  a family  of  geometrically  similar  reversers. 
In  passing  it  should  be  noted  that  the  error  made  in  the 
calculation  of  by  assuming  the  jet  leave  the  deflector  exactly, 
parallel  to  the  tar, gent  at  the  lip  of  the  deflector,  i.e,  by 
writing  /«  = 0,  as  ws  done  by  Schnurr  et  al,  would  be 
(assuming  Vc  remains  relatively  unaffected)  nearly  proportional 
to  ( I — cos  /I).  This  error  could  be  rather  large  particularly  for 
large  angles  of  deflection  and  its  correction  through  the  use  of 
another  arbitrary  parameter  such  as  the  “spillage  coefficient" 
introduced  by  Schnurr  et  al.':  is  not  justified  Evidently,  the 
correct  calculation  of  thrust  requires  the  determination  of  the 
deflection  angle. 

The  results  of  the  finite  clement  analysis  and  their  comparison 
with  those  obtained  experimentally  are  presented  in  Table  I 
together  with  those  obtained  through  the  use  of  the 
Belotserkovsky's  method.  The  experiment',  values  of  i]„  are 
somewhat  larger  than  those  obtained  numerically.  There  are 
several  reasons  for  this  difference;  the  two  most  important 
ones  being  the  Coanda  effect  (see.  e.g.  Ref.  35)  and  the  nozzle 
pressure  ratio.  The  Coanda  effect,  i.e,  the  tendency  of  a jet  to 
attach  to  adjacent  surfaces  (in  this  case  to  the  outer  surface 
of  the  nozzle)  because  of  the  entrainment  deprivation,  decreases 
the  deflection  angle  and  thereby  increases  the  reversed  thrust. 
In  'he  experiments  conducted  by  Steffen  et  al,34  the  outer 
surface  of  the  nozzle  \ as  streamlined  in  the  form  of  a boattail 
to  decrease  entrainment  and  thus  to  increase  In  fact,  as 
noted  by  Stcflen  et  al,  "the  pressure  reductions  on  the  boattail 
were  large  enough  to  account  for  as  much  as  2W/„  of  the 
. cvcrse-thrust  ratio."  Thus  the  relatively  large  differences  between 
the  ijk  values  experimentally  obtained  by  Steffen  et  al.  and 
those  predicted  numerically  are  primarily  attributable  to  the 
Coanda  effect. 

The  present  experiments  were  conducted  with  a straight, 
sharp-edged  nozzle  (r„  = 1.32  to  1.57  in,  Rc  = 2.5  in.)  and  the 
effect  of  the  jet  attachment  has  been  minimized.  This  has 
resulted  in  a closer  agreement  between  the  computed  and 
experimental  values  of  r/*.  In  passing  it  should  be  noted  that 
the  Coanda  effect  is  not  necessarily  desirable,  in  spite  of  its 
contribution  to  the  reverse  thrust,  for  it  may  cause  an  unstable 
llowfield  and  destructive  vibrations. 

The  nozzle-pressure  ratio  or  the  actual  velocity  in  the  nozzle- 
causes  variations  in  »/#  primarily  because  the  entrainment  reeds 
of  the  deflected  jet  and  hence  the  Coanda  effect  increase  with 
increasing  jet  velocity.  Consequently.  i;K  increases  (about  10 
percentage  points)  over  a range  of  nozzle-pressure  ratios  from 
1.7  to  3. 

Table  I shows  that  the  relative  spacing  of  the  nozzle  and  the 
reverser  size  significantly  affect  the  flow  reversal.  Evidently  ifK 
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increase  with  an  increase  in  hemisphere  diameter  and  reaches 
a nutximum  of  about  K(7’0  for  Ru  r„  = IS  and  sjr.  0 40 
This  nearly  corresponds  to  an  optimum  no//le-reverser  spacing. 
I \tensi\e  experimental  data  obtained  by  Steffen  et  al 34,J'’ 
with  various  types  of  reversers  clearly  show  that  at  spactngs 
greater  than  that  required  for  optimum  performance.  »/„  may- 
drop  as  much  as  20flo  Closer  spacings  do  not  noticeably 
affect  the  reverser  performance  but  result  in  decreased  flow 
rate  for  a given  total  nozzle  pressure  or  in  increased  pressure 
for  a given  How  rate  because  of  the  increased  blockage  or 
back  pressuring  effect  of  the  reverser 


Conclusions 

Though  the  equations  for  jets  impinging  on  curved  surfaces 
are  not  tntegrable  in  a closed  form,  the  behavior  of  their 
solutions  is  predictable  through  the  use  of  numerical  methods. 
I or  two-dimensional  cases,  the  Lcvi-Civita  method  provides  a 
variety  of  bucket  shapes  and  enables  one  to  round  off  the 
sharp  corners  of  buckets  otherwise  composed  of  straight 
segments.  Tor  these  cases,  the  finite  element  method  could  also 
provide  direct  solutions  for  prescribed  boundary  shapes. 

The  Belolscrkovsky  and  finite  element  methods  may  be  used 
with  confidence  for  the  analysis  of  an  asymmetric  I ai place 
field  where  one  or  more  parts  of  the  boundary  are  to  be 
determined  as  part  of  the  solution.  In  particular,  these  methods 
are  capable  of  predicting  the  idealized-performance  charac- 
teristics of  asymmetric  thrust  reversers  Evidently,  some  of  the 
practical  problems  aswcialcd  with  thrust  reversal,  such  as  the 
reattachment  of  the  jet  to  the  nacelle  of  the  engine,  hot-gas 
reingestion,  interaction  of  the  dellected  jet  with  the  ambient 
stream,  noise  intensification,  etc.,  require  additional  analytical 
and  experimental  investigation. 
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